Decay behavior of localized states at reconstructed armchair graphene edges 
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Density functional theory calculations are used to investigate the electronic structures of localized 
states at reconstructed armchair graphene edges. We consider graphene nanoribbons with two 
different edge types and obtain the energy band structures and charge densities of the edge states. 
By examining the imaginary part of the wavevector in the forbidden energy region, we reveal the 
decay behavior of the wavefunctions in graphene. The complex band structures of graphene in the 
armchair and zigzag directions are presented in both tight-binding and first-principles frameworks. 
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I. INTRODUCTION 

One-dimensional boundaries are a distinctive feature 
of finite-sized graphene, which means that an investiga- 
tion and understanding of the electronic properties of 
graphene edges are of particular importance. Even be- 
fore the successful production of graphene from graphite 
by exfoliation in 2004i, the edge state was shown to be 
a typical example of the manifestation of the topological 
characteristics of a bulk band^. The boundary or edge 
effects become more significant for practical applications 
as the size of graphene in a device becomes smaller. Real 
graphene edges are usually passivated with functional 
groups, transition metals, or hydrogen atoms depending 
on the specific purpose, but even without these exter- 
nal chemicals or elements, reconstruction of the graphene 
edge itself has been observed^. How deeply the edge- 
induced state penetrates the bulk determines the decay 
length and extent of the state localization. 

According to Bloch's theorem, the eigenstates of the 
single-electron Schrodinger equation in a crystal satisfy 
'^(r) = e'^^'^Uni^i'r), where link is a function that has the 
same periodicity as the crystal, n is the band index, and 
k is the wavevector. In the case of an infinite crystal, 
the Born-von Karman cyclic boundary conditions^ re- 
strict the wave vectors to real quantities. However, com- 
plex Bloch k vectors are allowed for finite crystals, and 
the complex band structure of a periodic system is the 
conventional band structure extended to complex Bloch 
k vectors. Near a crystal surface or interface, one can 
match a wavefunction with a complex k between the in- 
side and outside of the crystal region, and thus surface 
or interface evanescent states arise^— . 

The complex band structure concept can also be 
adopted for graphene edges; the properties of the edge 
states are closely related to the band structure of infinite 
graphene. The wavevector can be split into a component 
parallel to the edge, /c||, which is conserved during scat- 
tering, and a perpendicular component, k±. Then, for 




each real /c||, the dispersion relation E = E{k±) allows a 
complex k± = k' -\- in^ where k' and n are real. We refer 
to the imaginary part n as the decay parameter. 



FIG. 1. Model structures of reconstructed armchair edges: 
(a) pentagonal armchair edge with missing carbon atoms 
[ac(56)], and (b) reconstructed armchair edge in which two 
heptagons created by the Stone- Wales transformation share a 
side [ac(677)]. The bond lengths (in A) at the graphene edges 
are shown in the right panel. The gray and white spheres are 
carbon and hydrogen atoms, respectively. 



In this paper, we focus on reconstructed armchair 
graphene edges. A perfect armchair graphene edge has no 
localized edge state since its one-dimensional bulk Hamil- 
tonian Hk^^{k±) is topologically trivial for any k\\ valued. 
However, reconstructed armchair edges may have local- 
ized edge states because of the modification of the geome- 
tries and hopping properties. Using first-principles cal- 
culations, we show that the possible localized edge states 
have decay behaviors that are associated with the com- 
plex band structure of graphene in the bulk. In addition, 
we provide analytic solutions to the complex band struc- 
tures for both the armchair and zigzag directions using 
the nearest-neighbor tight-binding method. 



II. COMPUTATIONAL DETAILS 

We performed density functional theory calculations 
within the generalized gradient approximation (GGA) for 



the exchange- correlation functional. The Per dew-Bur ke- 
Ernzerhof (PBE) functional forrn^- was adopted for the 
GGA. Ionic potentials were described by the projector 
augmented wave (PAW)^^ method implemented in the 
Vienna Ab Initio Simulation Package (VASP)^^. Arm- 
chair graphene nanoribbons (GNRs) with ac(56)- and 
ac(677)-type edges^ on one side were calculated to mimic 
a semi-infinite graphene edge, and 24 and 12 k points 
were used to sample the Brillouin zone (BZ) in the edge 
direction in the respective edge geometries. A plane- wave 
energy cutoff of 400 eV was used for the structural re- 
laxation, which continued until the atomic forces were 
smaller than 20 meV/A. The size of the unit cell of the 
ac(677)-type edge was twice that of the ac(56)-type edge 
along the ribbon axis, as shown in Fig. [H The supercells 
contained 91 and 204 carbon atoms for the GNRs with 
the ac(56)- and ac(677)-type edges, respectively. 

To calculate the complex band structures in the primi- 
tive cell, we employed the Quantum Espresso packagei^. 
The Vanderbilt ultrasoft pseuopotentiaU^ was gener- 
ated through the Rappe-Rabe-Kaxiras-Joannopoulos 
(RRKJ)i^ pseudo-wavefunction construction scheme. 
The kinetic energy cutoff was 30 Ry. From the symmetry 
of the hexagonal cell, the net force exerted on each carbon 
atom was zero, and the lattice parameters were optimized 
by the Birch-Murnaghan equation of state^ii^. For the 
k point sampling, 24 x 24 x 1 k points were chosen in 
the Monkhorst-Pack schemeii. 

In the tight-binding approximation, the Hamiltonian 
of a graphene monolayer is expressed as 
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where XI ^^ .•> sums over only nearest-neighbor pairs, 

c\{ci) is the creation (annihilation) operator, and h.c. 
indicates the Hermitian conjugate. Here, t is the nearest- 
neighbor hopping energy (= —2.88 eV). 
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FIG. 2. (a) Band structure of a GNR with an ac(56) edge, and 
(b) electronic densities of two edge-related states. The labels 
A (B) and A^ (B^) indicate electronic densities in the same 
band at F and X, respectively, (c) Plane- averaged electron 
densities (red) and their semi- log plots (blue) of the GNR with 
the ac(56) edge showing an exponentially decaying character; 
the gradients of the fitted curves give the decay parameters 
in units of A~^. The position indicated on the horizontal axis 
is the distance from the hydrogen atom at the edge. 



III. RESULTS AND DISCUSSION 

The optimized structures of the two GNRs are shown 
in Fig. [H The ac(56) model shows a pentagonal re- 
construction of the armchair edge (with a connecting 
hexagon) that requires the diffusion of carbon atoms^. 
The ac(677) model^^— , however, shows a reconstruction 
in which two heptagons share a side (a carbon bond). Un- 
like the models of Koskinen et at? ^ our graphene edges 
are passivated with hydrogen atoms. 

Figure [2] shows the electronic band structure and 
charge densities at F and X in the armchair GNR with the 
pentagonal reconstruction, ac(56). Two bands labeled as 
AA' and BB' are seen in the forbidden energy band. Be- 
cause of the added carbon atoms, one additional band, 
labeled as AA', appears near the Fermi level. In terms 
of the tight-binding calculation, this corresponds to the 
inclusion of new basis functions. At the zone boundary 



(X), the edge states exist as deep levels in the forbidden 
energy region so that they are more localized than the 
edge states at the F point. The overlapping of a localized 
state with the bulk states in k-space is generally referred 
to as surface resonance^^; the localized state penetrates 
the bulk and couples strongly to the bulk states. Thus, 
the electronic states labeled A and B for the ac(56) model 
are extended GNR states with enhanced amplitudes near 
the edge. 

In our previous study^^ of the reconstructed armchair 
edge ac(56), we showed that the edge hopping energy 
to at the pentagon is approximately — 2 eV by using the 
maximally localized Wannier function method^l. In ideal 
graphene, the hopping energy t is —2.88 eV, which thus 
corresponds to boundary softening {to < t). Li et al^ 
showed that edge-hopping modulation may give rise to 
the edge state at the perfect armchair edge. 

We can estimate the decay lengths of edge states from 
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FIG. 3. (a) Band structure of a GNR with the ac(677) edge, 
and (b) electronic densities of two edge-related states. The 
labels C (D) and C^ (D^) indicate electronic densities in the 
same band at F and X, respectively, (c) Plane- averaged elec- 
tron densities (red) and their semi-log plots (blue) for the 
GNR with the ac(677) edge showing an exponential decay- 
ing character; the gradients of fitted curved give the decay 
parameters in units of A~^. The position indicated on the 
horizontal axis is the distance from the hydrogen atom at the 
edge. 



a calculation of the energy levels of a relatively narrow 
GNR because the interaction between the two edges of 
the nanoribbon quickly decreases as the width of the 
ribbon increases; the energies of the edge states then 
rapidly converge to the semi- infinite limit. The decay 
length Adecay(= 2ti/ n) is obtained from the decaying fac- 
tor e~'^^^, where r±_ is the perpendicular direction to the 
boundary. The estimated decay lengths are 1.23 and 1.59 
nm for states A' and B', respectively, as obtained from 
gradients of —0.510 and —0.394 A~^ in the semi-log plot 
of the decaying electronic densities in Fig. O The shorter 
decay length of state A' can also be seen in Fig. [2fb). 

Figure [3] shows the edge states near the Fermi level at 
the ac(677)-reconstructed edge. The reconstruction can 
be understood as an array of Stone- Wales defects, as 
this geometry is obtained by rotating particular carbon- 
carbon bonds by 90°. The Stone- Wales defect in car- 



bon nanotubes is known to be able to form two quasi- 
bound states^^ that are characterized as bonding and 
anti-bonding states in the rotated carbon dimer— . Fig- 
ure [31(a) shows that a nearly flat band, labeled as CC^ 
is unoccupied and may act as an acceptor level. In 
Fig. [31(b), the electronic charge densities of states C' and 
D have distinct anti-bonding and bonding characteristics 
in the rotated carbon dimer, respectively. Although the 
energy dispersion of the CC band is nearly flat, the decay 
lengths of the edge states vary significantly depending on 
k\\ in the edge direction. In fact, near the F point, the 
edge state is located inside an allowed energy band so 
that it is not a localized state due to mixing with the 
extended states. For state C^ the energy level is deep 
into the bulk energy gap. Consequently, its electronic 
charge density is strongly localized at the edge and has 
a decay length of 1.46 nm. In contrast, the energy level 
of state D' is so close to the valence band edge over the 
whole k\\ range that the decay length is relatively long 
(2.73 nm). The difference in the decay lengths can also 
be seen Figure [31(b). 

We can deduce the decay length of an edge state from 
its energy level when the Bloch wavevector k is assumed 
to be a complex number. The general procedure for con- 
sidering the complex wavevector is as follows: The crystal 
translational symmetry produces a wavevector k a good 
quantum number, and the bulk Hamiltonian can be de- 
coupled for each k. The secular equation 

det[E - i7(k)] = det[E - H(k\\ , A:^)] = (2) 

then gives the energy levels and Bloch wavefunctions of 
the crystal. In this case, k is regarded as a parame- 
ter. If we now also regard the energy and k\\ as param- 
eters and solve the secular equation with respect to /c^, 
then the secular equation becomes a polynomial of k^ 
in general^^. In line with the fundamental theorem of 
algebra, the polynomial has the same number of zeroes 
(including complex values) as the degree of the polyno- 
mial independent of the energy and k\\. This means that 
whenever there is a band edge where the real band begins 
to vanish, a corresponding complex band appears from 
the band edge. In view of the single-particle Schrodinger 
equation, the wavefunction at the boundary should be 
matched to a linear combination of bulk wavefunctions 
(including waves having complex wavevectors). When 
this linear combination does not include any wavefunc- 
tions with a purely real wavevector, it forms surface or 
interface evanescent states. 

Near the band edge, the behavior of the complex band 
is easily derived from an effective Hamiltonian. Usually, 
the real band has an approximately quadratic dispersion 
near its edge, and the effective Hamiltonian can be ex- 
pressed as 
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For a given k\\ = ko and E^ solving this Hamiltonian with 
respect to k± gives the complex band structure. When 
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4. Tight-binding complex band structures of graphene 
A:ii = 0, (b) k\\ = ^G\\, and 



in the armchair direction (a) k\\ = 0, (b) k\\ = jU\\, 



(c) k 



l^lh 



where G\\ 



^, and a is the carbon-carbon 



bond length (= 1.42 A), and the zigzag direction (d) k\\ =0, 

(e) k\\ = l^ii, and (f) k\\ = ^Gn, where G'n = 7=^- The 
hopping energy is set to —2.88 eV. 



E lies in the forbidden energy regime {E < 
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and the decaying parameter k± ex y^^^£^, where Eq = 

"2^ + £^0 is the band edge for a given k\\. 

The nearest-neighbor tight-binding model provides the 
essential features of the complex band structure and as- 
sists in understanding the overall structures of the com- 
plex bands; however, it does not provide an accurate de- 
cay length. We can solve the secular equation of a 4 x 4 
tight-binding matrix for k^ and obtain analytic solutions 
for specific ky and E values. In Fig. HI the tight-binding 
complex band structures are plotted for particular k\\ val- 
ues in both the armchair and zigzag directions. In both 
cases, the complex band has a quadratic shape near the X 
point (the band edge), the localization becomes stronger 
further away from the band edge, and reaches a maxi- 
mum deep inside the band gap. At any energy around 
the Fermi level, there are always four (two) 7r-bands in 
the graphene in the armchair (zigzag) directions. This 
feature is related to the different number of maximum 
scattering channels in the two directions because the en- 
ergy and crystal momentum along the direction parallel 
to the edge are conserved during electron scattering22. 
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FIG. 5. First-principles complex band structures of graphene 
in the armchair direction. Real and imaginary bands are plot- 
ted along the positive and negative axes, respectively. The 
lower hexagons are the BZ of graphene, and the yellow rect- 
angles denote the BZ if a doubled unit cell of graphene is 
considered for the armchair cell construction. The band struc- 
tures are plotted along line segments. From the left panel, (a) 



fc|l =0, (b) fell = iG||, 

and a is the carbon-carbon bond length. 



and (c) k 
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l^W^ 



where G^y = 3J, 



Complex band structures of graphene in the armchair 
direction calculated from first principles in a plane- wave 
basis are shown in Fig. [5] for three k\\ values {k\\ = 0, 
|G||, and \G\\)^ where Gy = |^, and a is the carbon- 
carbon bond length. In principle, any edge state in the 
bulk region can be represented by a linear combination 
of wavefunctions that correspond to a complex band of 
the same energy, and in practice, the dominant contri- 
bution comes from the first few bands with long decay 
lengths. Furthermore, in case of the graphene edge, the 
first complex band with the longest decay length is com- 
posed of TT-electrons and determines the decay length of 
the edge state. The energy levels of the edge states of 
the ac(56) edge at the X point (A' and B' in Fig. [2]) are 
—0.84 and —1.64 eV, respectively, and the corresponding 
decay parameters {k) are 0.255 and 0.207 A~^, as shown 
in Fig. [5fc). Because the charge density is given by the 
squared magnitude of the wavefunction, the decay pa- 
rameters {n) should be doubled for comparison with the 
gradients in Fig. [2l this gives decay lengths of 0.510 and 
0.394 A"^ for A' and B^ respectively. The small differ- 
ence is attributed to the finite size effect of the ribbon 
width. 



In the case of the ac(677) edge, both k\\ = \G\\ and 
k\\ = |G|| should be considered due to the doubling of 
the armchair unit cell. The energy levels of the C and D' 
states in Fig. [3] are 0.14 and —1.55 eV, respectively, and 
the corresponding decay parameters are 0.219 and 0.110 



A- 



1 at fc|| = iG| 



and 0.249 and 0.216 A'^ at fcii = hG\\ 
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FIG. 6. First-principles complex band structures of graphene 
in the zigzag direction. The yellow rectangle denotes the BZ of 
a doubled unit cell of graphene for the zigzag rectangular cell 
construction. From the left panel, (a) k\\ =0, (b) k\\ = ^G\\, 
and (c) k\\ = ^G\\, where G\\ = ^l~5 ^^^ ^ is the carbon- 
carbon bond length. 



Because the complex bands at /cy = jG\\ have smaller de- 
cay parameter, the decay lengths are determined by those 
bands. As mentioned earlier, the ac(677)-type edge has 
a lattice parameter that is twice that of the ac(56)-type 
edge, and the band edge corresponding to the ac(677) 
edge is at k\\ = \G\\. If the k values at k\\ = jG\\ is dou- 
bled, then agreement with the gradients of the charge 
densities (0.431 A-^ for C and 0.230 A-^ for D^ in 
Fig. ^c) is obtained. 

We would like to stress here that the relationship be- 
tween the complex band and the decay length is applica- 
ble to any graphene edge regardless of its chemical pas- 
sivation because only the energy dispersion of the edge 
state along k\\ is needed. In practice, since the calcu- 
lated energy quickly approaches the semi-infinite limit 
even for a relatively narrow ribbon, we can deduce the 
decay length of various graphene edge states from rela- 
tively narrow ribbon calculations. 

Grain boundaries in polycrystalline graphene, in a sim- 
ilar manner to reconstructed edges, also have topological 
defects such as pentagons and heptagons. At a grain 
boundary, at least two domains are matched, and the 
topological defects give rise to localized electronic states. 
In such cases, the decay lengths are affected by the crys- 
tallographic direction of the domain. The complex band 
structures in the zigzag direction for k\\ = 0, ^G|| and 
^G||, where G\\ = -?=-, shown in Fig. [6l have consid- 
erably different structures from those in the armchair 
direction (Fig. [5j). This unambiguously shows that the 
decay length of a localized state depends on the crys- 
tallographic direction of the domain when two distin- 



guishable domains are connected and a grain boundary 
is formed. We expect that the decay behavior of local- 
ized states originating from the topological defects at the 
grain boundaries is observable using scanning tunneling 
spectroscopy. 



IV. CONCLUSION 



In summary, we have shown that it is possible for local- 
ized states to appear in reconstructed armchair graphene 
edges using ab initio pseudopotential calculations. The 
edge state in the ac(677) model decays more rapidly than 
that in the ac(56) model at the X point. We have also 
presented complex band structures of graphene in the 
armchair and zigzag directions in both the tight-binding 
and first-principles frameworks. The extension of the 
conventional band structures to a complex band struc- 
ture provides information on the energy-dependent de- 
cay lengths of the graphene edge states. By comparing 
the shapes of the complex bands in the armchair and 
zigzag crystallographic directions, we revealed that the 
decay behaviors of the edge state are strongly related to 
the crystallographic directions. Our analysis indicates 
that our theoretical approach to understanding the edge 
states through the complex band structure is quite gen- 
eral and can be applied to any graphene-based structure 
with edges or grain boundaries. 



APPENDIX: TIGHT-BINDING MODEL FOR 

CALCULATING THE COMPLEX BAND 

STRUCTURE OF GRAPHENE 



If we construct the unit cell of graphene in the arm- 
chair direction as shown in Fig. [Tl^a), the Hamiltonian 
can be Fourier transformed, and each decoupled H{k) is 
represented as the following 4x4 matrix. 
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where k = {kx, ky)^Rx = (\/3a, 0)^ Ry = (0, 3a), and a is 
the carbon-carbon bond length. The order of the basis 
{i = 1, 2, 3, and 4) is shown in Fig. [Tja). If we define 
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secular equation for a given E and ky can be written as 
follows: 





FIG. 7. Doubled unit cell of graphene for the (a) armchair 
and (b) zigzag construction. The bases are indexed in the 
same order as the 4x4 matrix representation oi H{k) in the 
text. 
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(5) 



Because r 7^ 0, this equation is a fourth-order polyno- 
mial of r and has four solutions. Equation (5) is invari- 
ant under the operations r ^ ^ and ^ ^ r. If r^ is a 
solution to Eq. (5), then r*, — , and 4r are also solutions 
of the secular equation. If the solution has a nonzero n 
(= Im kx)^ then the four solutions correspond to left- and 
right-decaying states and their time-reversal pairs. The 
existence condition of the complex band can be derived 
if n is set to zero. 

To visualize the energy-dependent decay length, con- 
sider a special case of the band edge: ky = ^G\\ = f^. In 
this case, if we put s = — 1, the equation becomes 



1 o fE 

- + 3- - 
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= 0. 



(6) 



The solution of equation (6) is 
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and numerical solutions are plotted for the case of ky = 
0, g^, and ^ in Fig. Hl^a), (b), and (c), respectively. 

In the same manner, for graphene in the zigzag direc- 
tion, the Hamiltonian H{kx, ky) is given by 
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where k = {kx^ky), Rx = (3a, 0), and Ry = (0,\/3a). 
The order of the basis {i = 1, 2, 3, and 4) is shown in 



Fig. El^b). If we define 



^zk-R^ and s = e^^'^y , then 



secular equation for a given E and ky can be written as 
follows: 
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This equation has two solutions connected by an inverse 
relation. From this equation we can calculate the decay 
length of the edge states of the GNR. Because the edge 
state of a zigzag-edged GNR is a zero-energy mode, the 
/cjy-dependence of the decay length is obtained by solving 
the following equation. 



{2^2Re{s)} (r + 
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equation 
and 

are plotted for the cases oi ky = 0, 
Fig. m^d), (e) and (f), respectively. 
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